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Effects are investigated of CP symmetry violation in the decay of a scalar particle X (the 
Higgs boson) into two off-mass-shell Z bosons both decaying into a fermion-antifermion pair, 

—>■ Zi ZJ —>■ fifi / 2 / 2 - The most general form of the amplitude of the transition X —>■ ZiZJ, 
wherein the boson X may not have definite CP parity, is considered. Limits of applicability of the 
narrow-Z-width approximation used when obtaining differential widths of the decay under consid¬ 
eration are determined. Various observables connected with the structure of the amplitude of the 
decay X —>■ Z*ZJ are studied. These observables are analyzed in the Standard Model, as well as 
in models conceding indefinite CP parity of the Higgs boson. An experimental measurement at 
the LHC of angular and invariant mass distributions of the decay X —>■ ZiZJ —>■ / 1 / 1 / 2/2 can give 
information about the CP properties of the Higgs boson and its interaction with the Z boson. 
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I. INTRODUCTION 


a 


In 2012 the ATLAS and CMS collaborations detected [1| a neutral boson h with a mass of about 126 GeV. At 
the present time, detailed study of properties of this particle, called the Higgs boson, is an important task. The 
Standard Model (SM) Higgs boson is a state with = O'*"’*', and all the available experimental data about 
properties of the particle h are close to the corresponding theoretical predictions about the SM Higgs boson (see, 
for example, 2 t^)- In particular, the spin of the boson h is equal to zero or two, and many hypotheses in which 
the spin of h is two are excluded with probability 95% or higher [3|. At the same time, the situation may be more 
complicated. For example, some supersymmetric models predict existence of neutral bosons with negative or even 
indefinite CP parity 5l-l7|. 

The issue of the CP parity of the Higgs boson is also related to the search for CP symmetry breaking sources 
which are additional to the mechanism built into the Cabibbo-Kobayashi-Maskawa quark-mixing matrix. Such 
sources of CP violation, for example, in the Higgs sector, could help in explaining the known problem of the 
matter-antimatter asymmetry in the Universe [8|. 
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It has been suggested jsjllOl that the CP properties of the Higgs boson be studied by investigation of decays into 
two photons, h ^ via measurement of the polarization characteristics of thephotons. In Refs. [h| the decay 


to the photon and the Z boson, h — Z*^ —>■ //y, has been examined while 1^ [iSj study the decay to the photon 
and a lepton pair, h —^ . In these papers it has been shown that the “forward-backward” escape asymmetry 

for the final fermions carries information about the CP properties of the h boson and physics beyond the SM. 

Investigation into the decay of the Higgs boson into two Z bosons with their consequent decay to fermions is 
another opportunity to ascertain the CP properties of h. Such a cascade decay wherein the final fermions are 
leptons, along with the two-photon decay channel, has allowed the determination [lj| of the mass of the particle h 


with the highest accuracy. In Refs. 


14 


17 theoretical distributions of the decay h —>■ —>■ /i / i /2 f 2 have been 


studied at various values of the spin of h and in case of various CP properties of this boson. In 


14| it has been 


reported what properties of experimental distributions testify about a particular spin and a particular CP parity 
of h. In 1 ?tI 417I| asymmetries measurement of which allows clarification of the mentioned properties of the Higgs 


boson are suggested and investigated. Finally, papers 


18l | put forward various methodologies on getting constraints 


on the Higgs boson couplings to ZZ, W~W^, 77 and Z 7 from experimental data. 

Besides, various theories with spontaneous breaking of the conformal invariance (for example, theories of tech¬ 
nicolor) assume the existence of one more neutral zero-spin particle which interacts with the gauge bosons - the 


m 


2(h 22I| it has been shown that the variant in 


dilaton. At present, the mass of the dilaton is not determined, but according to estimates performed in Ref. 
in some models the mass can exceed lO'*’ GeV. Along with that, in 
which the boson h is the dilaton is not excluded. 

In order to clarify the CP properties of the particle h and the hypothetical dilaton we consider a neutral particle 
X with zero spin and arbitrary CP parity. We examine the decay X —>■ Z^Z^ —>■ /1/1/2/2 in case of the non¬ 
identical fermions, /i fy fy, and study in detail the differential width of this decay with respect to the three angles of 
the fermions in the helicity frame and with respect to the invariant masses of the fermion pairs /i/i and fyfy- The 
most general X —> Z* vertex, which generalizes the corresponding SM vertex and contains a term corresponding 
to the negative CP parity of the particle X, is used. 

We also find limits of applicability of the narrow-width approximation for the Z boson for the presented calcu¬ 
lation of differential widths of the given decay. By means of this approximation we derive a formula for the total 
width of the decay X ^ Z\Z^ — >■ /1/1/2/2 (the formula is valid also in case fy = f 2 ) and a formula for the total 


width of the decay h ^ ZfZ^- These formulas are more general and more precise than those obtained in Ref. 231. 

Next we find observables connected with the structure of the amplitude of the decay X —>■ Z^Z^- The formula 
for the fully differential decay width contains nine coefficients related to the amplitude X —)► ZIZ^. For each of 
them one or two observables linear in this coefficient are dehned. Note that some of these observables, as well as 


different ones, have been studied in 


15-17 


24| , however we also obtain new experimentally measurable quantities 
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and analyze the dependences of the observables on the mass of one of the Z bosons (Z|) in much more detail than 
it has been done in the mentioned papers. This analysis is carried out within the framework of the SM as well as 
in certain SM extensions wherein the boson h is a mixture of a CP-even state and a CP-odd one. Measurement 
of the suggested observables at the LHC can yield important information about the CP properties of the Higgs 
boson and its interaction with the Z boson. 


II. FORMALISM FOR THE DECAYS A ZlZl / 1 / 1 / 2/2 

A. The amplitude of the decay X ^ ZIZ 2 and the fully differential decay width for X ^ ZIZ 2 / 1 / 1 / 2/2 

Let us consider the decay of a neutral spin-zero particle X with arbitrary CP parity into two off-mass-shell Z 
bosons (Z^ and Z|) each of which decays to a fermion-antifermion pair, /i/i and / 2 / 2 , 


X —>■ Z1Z2 —>■ /1/1/2/2, 


( 1 ) 


where mx > 2(to/j + 2 ) (to satisfy the law of conservation of energy in a rest frame of X), mx is the 
mass of the particle X, mf. is the mass of the fermion fj. We will consider this decay at tree level. If 
mx G (477i;,,2mt] {mb is the mass of the b quark, mt is the mass of the t q uark), which holds true if A = h. 


then fj = e , /i , r s, b. If mx > 4mt, which is possible 

the top quark as well. 


19j if X is the dilaton, then fj can be 


From the energy-momentum conservation we find that ai and 02 (oj is the mass squared of the boson Z*, i.e. 


the invariant mass squared of the pair fjfj) lie within limits 


4ml < oi < {mx - 2m 4to1 < 02 < {mx - v^)^- 


( 2 ) 


The amplitude Ax^z*z* (Ai, A 2 ) of the decay of X into Z* and is equal to [I5l - ll7l . l24| 

2lx^z*z*(Ai, A 2 ) =2\fj^^m%(az{e\ ■ e*^) P ^^{el ■ {pi +P2)){e*2 ■ {pi +^ 2 ))+ 

V m^ 

;W(K+P^)(K-pD(e?r(e2)*), 


cz 


(3) 


where Xj, e^, pj are respectively the helicity, the polarization 4-vector and the 4-momentum of the boson Z*, Gp 
is the Fermi constant, mz is the mass of the Z boson, az, bz, cz are complex-valued dimensionless functions of oi 
and 02 , is the Levi-Civita symbol (£0123 = !)■ Note that at tree level 

- if A is the SM Higgs boson, then az = 1, bz = cz = 0; 

if the CP parity of A is -1, then az = bz = 0 and cz 0; 


if the CP parity of A is indefinite, then 7 ^ 0, 7 ^ 0 and/or bz 0, cz 0. 
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Calculating the Lorentz-invariant amplitude (Ai, A 2 ) in a reference frame in which Pi + P 2 = 0, we 

derive that 

X^mx,ai,a2)\ 
m\ h 


(- 1 ,- 1 ) = 2 y V2GFm% i az — cz 


A5(TO|-,ai,a2)\ 


(1,1) = 2y V2GFm% ^z + cz 

(Ai, A2) = 0 , Ai ^ A2, 


m 


X 


(4) 


where the function X{x, y, z) is defined in the standard way: A(a;, y, z) = — 2xy — 2xz — 2yz. 



Fig. 1: The kinematics of the decay X ^ Zi —>■ fifif 2 f 2 - The momenta of Zi and Z 2 are shown in a rest frame of X, 

the momenta of /i and /i (/2 and / 2 ) are shown in a rest frame of Zl (Z^). 


To describe the decay (HD, let us introduce the following angles (see Fig. [Ij: 0i (O 2 ) is the angle between the 
momentum of Z* (Z 2 ) in a rest frame of X and the momentum of /i (/ 2 ) in a rest frame of Z* (Z 2 ) and (p is 
the azimuthal angle between the planes of the decays Z^ fifi and Z 2 / 2 / 2 - Further we go into the case of 
the non-identical fermions, /i / 2 . Using the helicity formalism (see, for example, 25|), we obtain that in the 
approximation of the massless fermions, rrif-^ = nr /2 =0, the differential decay width of ([I} with respect to oi, 02 , 
01 , 02 , f appears as follows: 


V2Glm^ 


F"‘’Z (2 i„2\/2 2 N “ 2)0102 


daida2d9id02d(p (47r)®TO|. ^ D{ai)D{a2) 

X sin0i sin02[(|A|| + |Aj_|^) ((1 -|- cos^ ^i)(l + cos^ 82 ) + AAf^Af^ cos0i cos 02 ) 

-I- 4|24o|^ sin^ 0i sin^ 02 — 4Re(A|2lj_)(24/j cos0i(l -h cos^ 62 ) + Af^ cos02(1 -|- cos^ 0i)) 
-I- 4-\/2sin0i sin02((Re(AgA||) cosy? — Im(24QAj_) sinyj)(-|- cos0i cos02) 


— (Re(24QAj_) cos 1/5 — Im(24g24||) sinyj)(T/j cos02 -I- Af^ cos0i)) 

-I- sin^ 01 sin^ 02((|A|| p — |A_lP) cos2i/j — 2 Im(AJ|24j_) sin 2:/))], (5) 


where o/ is the projection of the weak isospin of a fermion /, Vf = aj — 2^ sin^ 0w, qj is the electric charge of the 
fermion /, e is the electric charge of the positron, 9w is the weak mixing angle, D(ai^ 2 ) = (“ 1,2 ~ Fn\'f‘ + (rnzl'z)'^, 
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is the total width of the Z boson, Ar = . 

^ ’ J ai+vi ’ 


^ 11 ( 01 , 02 ) = 


^±( 01 , 02 ) = 


Ao(«li 02 ) = 


Ax->-z*z^(i,i) + Ax->.z*z*(-l, -1) _ 

23 y/Cpm"^ 

Ax^z*z*{i, 1 ) - Ax-^z•z•{-^, - 1 ) 

2iy/GF'm^z 

Ax^z*z*{0,0) ^ — ai — 02 

2i G pTn^z ^ ‘2,yJa\a-2, 


V2az, 


/- A^to^, 01 , 02 ) 

y zcz - o - 

mj^ 

A(TO^,ai,a 2 ) \ 

^ 4m|-^aia2 J ' 


( 6 ) 


Futher the approximation ruf^ = ^/2 = 0 is used. Using Eq. @, one can connect the ratios of quantities 
l^oPj 1^111^ + l^-Lp) 1^111^ ~ l^-Lpi R-g(^o^II ), Im(^Q^||), Re(AlQAj_), Im(^Q24j_), Re(24y24j_), Im(Ay Aj_) to + 
1^11 P + \A±\'^ with functions of oi, 02 which can be measured in experiment. We will call these ratios the helicity 
coefficients of the decay X —>■ Z^Z^- 


B. A differential width . F 

aCLl ClCL ‘2 


The number of the decays 


^ 1 *^ 2 * 


^l "^2 ^2 — 6 )/^); 


(7) 


a 


detected in the ATLAS experiment (2| wherein the invariant mass of the four leptons was in the interval [120 GeV, 
130 GeV], is equal to 32. The number of the decays ([71) detected in the GMS experiment |3j in which the four-lepton 
invariant mass was within [121.5 GeV, 130.5 GeV], is equal to 25. In view of the insignificant amount of data, at 
the present time an experimental dependence of the distribution ^ daidatdeid 92 dip (7" is the total width of the decay 
O) for any of the decays © is not available. Let us consider differential decay widths of 0 with respect to four 


and fewer variables. Integrating Eq. ([51) with respect to 0i, 02 , we obtain 


(PT 


V2G' 


daida2 9(27r)^m^ 




A^m|., 01,02)0102 


D( 0 i)I?( 02 ) 


E 


( 8 ) 


P= 0 .||.-L 

It follows from Eqs. ([5]), (© that the dependence of the differential width on az, bz, cz boils down only to 

the dependence on \az\, |fcz|, \cz\ and on cos(arg 62 — argoz). 

The available experimental data on properties of theparticle h are close to the corresponding theoretical predic¬ 
tions about the SM Higgs boson (see, for example, 2 t 4)- That is why ahz ~ 1, bhz ~ 0, Chz ~ 0, where 


ahz = az\x=h, bhz = bz\x=h, cuz = cz\x=h- 

In Fig. H we show the differential decay width ([ 8 ]) for A —— >■ l^l 2 I 2 {Ij = e,/i,T, /i ^ I 2 ) as a 
function of yja-i in the SM for |az| = l, 62 =c^ = 0 and mx = w/j, where rrih is the mass of the Higgs 
boson h. The range of ^/a 2 hr this plot is determined by the inequalities ([2|) in the approximation of the 

massless fermions. In calculations and when plotting graphs the experimental data listed in Table [Tj are used, and 
sin^ 9w = 1 — 'fn\rl'<dn?z^ where mw is the mass of the W boson. 
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10 

d^T .,,1 

' -,10 -r; fel = l,*z = 0,C2 = 0 

dij di2 G«V-* 

0 

100 

G«V 

0 

Fig. 2: The dependence of the differential decay width of X —>■ —>■ {Ij = e,/i,r, Zi 7 ^ I 2 ) on and 

in the SM for mx = rn^,. 

Table I: Values of the Fermi constant, of the masses of h, Z, W and of the total width of Z 

Gf = 1.1663787(6) x 10 "® GeV^ 
mh = 125.7(4) GeV 
mz = 91.1876(21) GeV 
mw = 80.385(15) GeV 
Tz = 2.4952(23) GeV 

As one can see from Fig. [2j in the SM the function has peaks at y^ = mz and y/o^ = mz^ resulting 

from the quantities D{ai) and £>( 02 ) in ([5]). 

Let us calculate the ratio of a typical value of in the SM on the peaks to its typical value in an area in 

which y/oi and y/^ significantly differ from mz (we will call this area “plateau”). As indicative values of ^/oi and 
y/oi’ on the peaks we take y/oT = mz, y/a 2 = — rnz) and y/oi = \{mh — mz), y/o^ = mz (see ([ 2 ])), and values 

on the “plateau” are chosen y/oi = y/^ = \mz. It follows from ([5]) that in the SM for any /i, /2 values of 
at y/oT = mz or y/o^ = mz are approximately 100 times as great as values of this function on the “plateau”. 

If mx 7 ^ mh but just greater than mz, then y/oi and/or y/a 2 can be equal to mz (according to (H])), and, 

t2 

consequently, in this case the behavior of the function in the SM is similar to that in case mx = mh- That 

is why for any mx > mz and for any final fermions the differential width in the SM has a sharp maximum 

at yjoi = mz or y/o^ = mz- Therefore, if jazj Ri 1, 6 ^ ~ 0, Cz Ri 0 (which is the case of a small distinction between 
the couplings and their SM values), also has a sharp maximum at y/^ = mz or y/^ = mz, provided that 





mx > mz- 
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C. Limits of applicability of the narrow-Z-width approximation 


In Refs. 


27 


29| the accuracy of the narrow-width approximation has been studied for calculation of the total 


widths of various decays along with the total and differential cross sections of various processes. It is shown that in 
many cases (especially for processes beyond the SM) this approximation is not applicable. In this connection the 
question arises whether the narrow-Z-width approximation is applicable for obtaining the differential width by 

l2 

means of integrating In this subsection we hnd the interval of all the a 2 -values for which the approximate 

integration is valid. 

We consider the tojv- values such that mx > rnz and the dependences of 02 ( 01 , 02 ), 1 ) 2 ( 01 , 02 ), 02 ( 01 , 02 ) such 
that for any fi and /2 has a sharp maximum when = m 2 or y^a 2 = m 2 (an example of such dependences 

is |a 2 | ~ 1, 62 ~ 0, cz ~ 0). Then while calculating the differential width ^ one may use the narrow-Z-width 
approximation: 


dr 

da2 




da I 


d^r 

daida2 


(mx--v/a 2 )^ 


doi 


m 2 r 2 


d(oi - m|)/(oi, 02 ) 


m 2 r 2 


/(m|, 02 ) V S (0, mx - m 2 - A], 


(9) 


where A is some positive quantity and 

- E ( 10 ) 

P=0.||,-L 

A > 0 since in Eq. one may use the approximation d(oi — m^)f{ai, 02 ) only when ^^02 < mx — 

i2 'pi 

m 2 , because if approaches mx — mz, the peak of at y^oT = mz gets less sharp and at y/d 2 = mx — mz 

the peak disappears (see Fig. [5] and Eq. ©)■ However, the derivation Q does not allow one to estimate the 
accuracy of the formula ^ Ri mzTz /(^z W 2 ) at a given value of ^/a 2 ^ and for this reason it is not clear what 
value of A should be chosen. 

To clarify this point, let us derive the formula for ^ in the following way: 


/(ai,a 2 ) = 




mV 


9(27r) 


(“.tl + + ^f2) 


2 . A 2 (m|-, 01 , 02)01 




D{a2) 


dr 

do 2 


{mx—x'f^Y m2+£2 

f ^ d^r f , d^r 

I ddi 


daida2 


0 

m|+e 2 


/ 


doi 


doido2 


/ 


da-I 


f{ai,a2) 


(oi - m|)2 -t {mzTzy 


da\ 


/(m|, 02 ) 


(oi - m|)2 (m2r2)^ 


arctan —%—h arctan 

- 02 ) 

m 2 l 2 


-/(TO 2 W 2 ), 


( 11 ) 


m 2 r 2 ‘ 

where ei and £2 are some positive quantities such that m 2 r 2 “C Sj <C m^, the variable 02 takes values in the 


interval f 0 , (jnx — \/m'z + £ 2 ^ 






























One of the approximations used in Eq. m is the switch from the integration over an interval ( 0 , {mx — 
to the integration over an interval (m^ —ei, +£ 2 )- Thus, m\ — £i has to be greater than or equal to 4m(which 
holds true since ei <C m^) and m^ +£2 has to be less than or equal to {mx — i-e. 02 < {mx — -^/ml + £ 2 ^ . 

The latter inequality restricts the interval of all the a 2 -values for which these approximations are applicable. 
Consequently, in order to apply them for as long an interval of 02 -values as possible, one should use the minimal 
£ 2 -value at which the approximations are valid. 

While obtaining (fTTll . we also used an approximation A k, x {A = arctan + arctan ). Let us define 

2 _ 

£1 as £1 = mzy/mz^z (so as ^^p^ = The values of quantities A and mt — \J m^ + £2 which are listed in 

Table [IT] specify for the considered £ 2 -values the accuracy of the approximation A Ki tt and the maximal value of 
at which the narrow-Z-width approximation is applicable in case X = h. 


Table II: Values of A (ei = mzy/mz^z) and of mj, — y/rri^ + £2 at various values of £ 2 - 


£2 A mh — y/rn% + £2 (GeV) 


0 O.IStt 

34.51 

mzTz O.TOtt 

33.27 

2mzVz O.SOtt 

32.05 

2>mzVz O.SStt 

30.84 

AmzVz O.STtt 

29.65 


According to Table HIl if £2 < 3mz^z, then A < O.SStt and, in view of the big difference between A and tt, we 
will not apply the approximations (fTTl) for such values of £ 2 . Hence we will use £2 = imz^z- It follows from (fTTl) 
that 


dr 

da2 




m'n 


9 ■ z 








A 2 (m 


m 


z^ 


02)02 


D{a2) 


P=0,||,-L 


|24;i 


V y/a2 G 


^ 0 , mx 



( 12 ) 


where A'^ = Ap{m%, (p = 0 , ||, T). 

Note that in Refs. Ij,ll6|,[l7| when plotting dependences of on y^, formulas for which correspond to 
(fT^ have been used, but these graphs have been plotted for y/a 2 < mx — mz, despite the fact that Eq. (ITT]) is not 
valid at £2 = 0 (see Table HI!) . and, therefore, the plotted dependences significantly differ from the true ones in the 


interval y^ G {mx — \J 2 ,mx — mz]- 
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D. An inequality constraining c'y^z from CMS data 


a 


According to |3|, 

(j{pp h) nh^zlz-^ii) 

cTSMipp ^ - 

where cr(pp —^ h) is the cross section for production of h in pp collisions, 


TsM{h-^ziz‘^Ai) ~ 0-93to:23(stat);o Jg(syst), 


(13) 


r(h ^ ZIZ^ 4/) = T{h ZIZ^ 4e) + r(h ^ ZIZ^ 4/x) + T{h -)■ ZIZ^ 2e2/i) 
= 2T{h^ ZIZ^ 4e) + r(h ^ ZIZ^ 2e2p), 


(14) 


F/j is the total width of the boson h, asuipp —>■ h), TsM{h Z^Z^ 41), ThSM are the predictions of the SM for 
respectively (T{pp —>■ h), r{h —>■ Z^Z^ —>■ 41), F/j at rrih = 125.6 GeV. Obtaining (fT51l . the CMS collaboration has 
combined data from pp collisions corresponding to an integrated luminosity of 5.1 fb“^ at a center-of-mass energy 
y/s = 7TeY and 19.7 fb"! at ^ = 8 TeV. 

We consider the case in which the functions \b'i^z\> Whz\^ cos(argl )(^2 ~ arga'^^) do not depend on 02 . Here 

we define 


O'hz ='^hzim^, a2), b'l^z = bhz{m^, a2), c'j^z = Chz{m^,a2)- 


Then using the approximation 


a{pp -)> h) (Jsm{pp h) 


F?i 7'hSM 

and Eqs. (ITSl) (within one standard deviation), (lAlOII . (IA13I) (see Appendix lAl) . we derive the relation 
|a;^|2 +0.015 |5',^|2 + 0.177Re(aCAz) +0.037e [0.68, 1.22]. 


(15) 


(16) 


While obtaining (1161) we plugged the central values of rrih, mz, F^ listed in Table U into Eq. (lAlOE Note that 
the latter equation is derived at tree level and without allowance for the interference term connected with the 
permutation of the identical fermions in case /i = / 2 . The interference contribution to T{h —?► Z*Z 2 4Z) at tree 


level is expected to be negligible since in the SM at nih = 140 GeV it amounts to 2.99% (see Table 1 in Ref. 30j). 
Using the data of Table iHll and considering two sigma errors where available, we obtain that at y/s = 8 TeV 

a{pp h)/Th 


S (0.17, 00 ), 


(17) 


O'SAliPP —>■ h)/ThSM 

which means that the approximation (11511 does not contradict the experimental limits. 

Moreover, assuming that all the couplings of the Higgs boson except for ahz, b^z and Chz are equal to their 
SM values, we can verify m- In this case the only anomalous contribution to F/i comes from F(/i —>■ Z*ZJ), 


which makes up, in the SM, only about 2.81% 


34| of the total Higgs boson width, and therefore F^ is unlikely to 
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Table III: Experimental and theoretical results for the total production cross-section of the Higgs boson in pp collisions and 
for its total width. 

<j{pp —> ft) = 33.0 ± 5.3(stat) ± 1.6(syst) pb at ^/s = 8 TeV [3^ 

(^sm{pp h) — 22.09 pb (uncertainties not available) at ^/s = 8 TeV 
Tj, < 22 MeV at 95% confidence level (CL) [3^ 

ThSM = 4.15 ±0.16 MeV |34l 


substantially differ from its SM prediction. Besides, the inequality (fTBl) means that |r(ft Z^Z^) — ^SM{h 
ZtZ 2 )\/TsMih —>■ ZIZ 2 ) G [0,0.32] because its left-hand side is 

r(ft ^ z*z* 41) _ r(k z*z*) 

rsM(^ z*Z2 —>• 41) rsM(^ ■^ 1 ^ 2 ) 

(see (fT^ . (lAlOII ). For this reason (fTHl) implies that the relative change of is less than 2.81% • 0.32 Ri 0.90%, and, 
consequently, ( 1161 ) is consistent with the approximation F^i Ri Fji sm • 

The dominant contribution to the Higgs boson production cross section <jsm(j>P h) comes from the gluon 
fusion process gg —^ ft, which is independent of the hZZ vertex. The processes involving the hZZ interaction, i.e. 
the Higgs-strahlung Zh and the Z boson fusion, constitute much less parts of asmipp ft)- Specifically, at ^/s = 


8 TeV they can be estimated as 0.41 pb and 0.70 pb respectively [S^- The total production cross section at this 
energy is 22.09 pb (see Table 11111) . so the processes of interest contribute about 5% of the total cross section. That 
is why it seems improbable that the couplings ahz, b^z and Chz provide a significant difference between a{pp —>■ ft) 
and <jsm{j>P ft)- However, a derivation of the dependence of the total production cross section on the hZZ 
couplings would require a separate study. 

Summarizing the discussion of the approximation (1151) . we can infer, firstly, that it is consistent with the available 
data 3l|, 1 3 31 ] and, secondly, under the assumption that the only anomalous Higgs boson couplings are related to 
the hZZ vertex, Eq. (fT51) is most likely to be valid due to the small contributions of the hZZ vertex to F^, and 
cr{pp —>■ ft). 


E. Constraints on c'y^z 

The inequality (1161) constrains the whole six-dimensional space formed by the real and imaginary parts of the 
couplings ftl v a nd to the set of ellipsoids allowed by (ITSl) . Note that a similar interpretation has been 
suggested in Ref. [35 1. 

From (iroi) it follows that the variant = 0, = 1 (negative CP parity of the boson ft) is excluded. 

Now let us find constraints on the values of assuming that is taken from the SM, i.e. ~ ^ 
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or a '^2 = 1- Then 


Whz\ — 1 — 0 Whz\ ^ [0i 2.44]; 

a'f^z = 1 c'f^z = 0 lm6^2 = 0 ^ 6'^^ G ([—12.66, —9.31] U [—2.22, 1.14]), 


a 


4^ = 1 and = 0 and = 0 ^ [—3.84, 3.84]. 


(18a) 

(18b) 

(18c) 


Let us compare (fT51) with the hZZ coupling constraints obtained by the CMS 3g and ATLAS [37| collaborations. 
For this purpose we first express our XZZ couplings in terms of the CMS ones Oi, a^, 0,3 (we denote oi, 02,03 from 
361 ] as hi, 0 . 2 , 03 to avoid confusion): 


az = a ( 3,1 - exp(i(/>Ai ) 


Ui + 02 — 01 — 02 

Af m| 


02 


bz = - 2 a—^ 


02 , 


Cz = -la- 


''X 


0-3, 


(19a) 

(19b) 

(19c) 


where a = aov/2, ao is the proportionality factor of the amplitude A{HZZ) of the transition X ^ Z^Z^ (see 
Eq. (1) in [361), V = Ilyjy/2GF is the vacuum expectation value of the Higgs field, Ai is a scale of physics beyond 


36| 


the SM, (pAi is the phase in the term with Ai. In general hi, 02 , 03 may depend on oi and 02 , however in 
they are set to be constant. The ATLAS XZZ couplings a, ksm, khzz, i^azz are related to the CMS ones in the 
following way: 

a ^hi - exp(i^Ai)—= KSMCOsa, aa 2 =-^khzz cos a, aa 3 =-^kazz sin a, ( 20 ) 

where A is the EFT energy scale. Note that comparing the Lagrangian (1) in [37| with the one describing the 


interaction of the SM Higgs field with ZZ and W~W~^, one can deduce that the coupling gnzz from (1) in 
equal to 2m\jv. In 37| the couplings a, ksm, khzz, kazz are considered constant and real. 

In Refs. 3al37j 95% CL allowed regions for hZZ couplings are reported (see Table HV)l . Note that 


33 is 


Khzz = -^khzz, 
Khzz 


kazz = -^kazz, 


m% bz 


kazz 


tana = 


2 im% Cz 


KSM 2 m\az + {'m\-oi- 02 )bz' ksm 
and in the limit Ai —>■ 00 the CMS and ATLAS ratios coincide: 


2 m\ az + {m\ — ai — 02 ) bz 


( 21 ) 

( 22 ) 


02 .. Khzz 

— = lim -, 

oi Ai-^-oo Ksm 


03 


= lim 


kazz 


oi Ai-^-oo y Ksm 


tana 


(23) 


Following 37|, we assume the ATLAS hZZ couplings to be constant. Then considering the case khzz = 0, we 
find that our couplings Ohz, bhz, Chz are constant as well (see (fnil) . (EUl) '). and using (I18al) we obtain an allowed 
interval for kazz ta.n a/ ksm (see Table IVTl . However, in case khzz = 0 the results (Il 8 bll and (I18cl) only show that 
h may be the SM Higgs boson, and thus they do not constrain any hZZ couplings. 
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Table IV: The CMS 


36 1 and ATLAS 371 95% CL allowed regions for hZZ couplings. 


The last row shows the conditions 


under which these regions have been derived. 



CMS 


ATLAS 


Q3 

KHZZ -KAZZ tana 

ai 

ai 

«5M «5M 

[-2.28,-1.88] U [-0.69, 00 ) 

[-2.05, 2.19] 

(-0.75, 2.45) (-2.85, 0.95) 

Im^ = 0, dAi = 0 or TT ImS2. 

ai ’ ai 

= 0, (^Ai = 0 or TT 

KAZZ = 0 Khzz = 0 


If ^Hzz 7 ^ 0, then a^z acquires a dependence on the invariant masses squared oi and 02 , and therefore the 
constraints m and (1181) get invalid since they have been derived under the assumption that \a'i^z\ I® independent 
of 02 - Therefore to constrain the ATLAS couplings in case khzz 7 ^ 0, we start with Eqs. ([Till) and (fTKl) . which 
demonstrate that within one standard deviation 


r{h -7 ziz^ -7 4 /)/rsM(/i -t ziz^ -7 4/) e [0.68,1.22]. 


To obtain T we have to calculate the integral (EH) for depending on 02 - Taking into account the limits of the 
integration, we substitute 02 with (rax — rrizY in the expression for a'^ (see (I19al) . (I^Hll l and therefore derive 
Eq. (lAlOl) where az has the expression (I19al) with 02 = (mx — rnz'f' It means that if khzz is not zero, we 
may use (ITHl) and (Il8bl) . (I18cll with determined by Eq. (I19al) where 02 is replaced by (mh — mzY This 
conclusion allows us to constrain khzz/ksm and ImkHzz/^s ksm, as one can see in Table fVl 


Table V: Our allowed regions for the ATLAS hZZ couplings. The last two rows show the conditions under which these 


regions have been derived. 


Khzz 

K-SM 


tan a 

KSM 


Im Khzz 

Re KgM 


[-2.38,-1.89] U [-0.24,1.13] 

02 = (mh — mzf' 1^ in (I19al) . 
o-'hz = «Azz sin a = 0, hnKHZz = 0 


[-1.28, 1.28] 

jo^zi = 1, rLHZz = 0, Im- 


•^AZZ _ g 
KSM 


[- 1 . 01 , 1 . 01 ] 

02 = (mh — mz)^ in lll9al) . 

(ihz = b kazz sin a = 0, ^bkhzz = 0 


Note that the results (ITB)) . (01^ along with the regions shown in Table IVl are estimated with consideration of 
the one sigma interval in (1131) . with the approximation (1151) . the central values of ruh, rnz, Tz from Table |T] and 




36 

37 


Refs. 


37 [and our ones. In addition, we present an allowed interval for the ratio Ivcikuzz k,sm unconstrained in 


Ml37|. 
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Fig. 3: The distribution as a function of for the decay h ^ ZIZ 2 —> / 1 / 1 / 2/2 in case \a'i^z\ = li ^!iz = 0) = 0 

(solid line); = 1, = 0i |c/izl = 0-5 (dashed line); = 1, h'f^z = —0.5, = 0 (dash-dotted line). 

We choose the following sets of values of a'^^, and c'^^: 

Whz\ = 1) ^'hZ = Oj c!hZ = Oj 

^'hZ = 1 ; ^'hZ = 0 ; ^'hZ = 

^'hZ = 1 ; ^'hZ = 0 , = ^. 5 *, 

= —0.5, Cj^z = 0 (24) 

and 

oJtZ = li ^'hz = —0.5f, = 0, (25) 

which are consistent with the constraints m- The sets (HI and (1^51) will be used for examination of further 
results. 

Regarding the selected values in (12^ and (1^51) one should mention that even in the SM the couplings b^z and Chz 
acquire small values due to electroweak radiative corrections where Imbhz and Imchz come from the absorptive 
parts of the corresponding loop diagrams. In Eqs. (P^ . (1^51) we assume that the hZZ vertex may be significantly 
modified by physics beyond the SM. 

It is of interest to study the distribution ^ function of y/a 2 for various sets of a'^, b'^, c'^. Here 

a'z = az{m%,a 2 ),b'^ = bz{rn‘^zja 2 ),c'^ = cz{rn'^,a 2 )- In accordance with ([8|), the function is independent of 

the final fermion state. Figure |3] shows this observable in case X = h. 

As one can see from Fig. |3l the function is sensitive to and almost insensitive to For this reason, 
having measured this distribution with sufficient accuracy, one can get significant constraints on the values of b'f^^. 
However, one should keep in mind that this conclusion is obtained for the case in which Whz\ 
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cos(arg 6'^2 ~ independent of 02 , and their a 2 -dependence can considerably modify the dependence of 

In Sec. Ill FI we develop methods of getting constraints on the dependences of a'^, b'^, c'^ on 02 - 


F. Connection between the helicity coefficients of the decay X ^ ZiZ^ and observables 

Let us consider now arbitrary dependences of 02 ( 01 , 02 ), & 2 (ai, 02 ), 02 ( 01 , 02 ) such that the differential width 
daida 2 deid 02 dip ^ sharp maximum as a function of oi and 02 at y^oi = mz or ^^02 = for any /i, / 2 . 
From Eq. dSD, using approximations analogous to those used when deriving the formulas (mD, IMl), we carry out 
integration over oi and some of the angular variables. Then we obtain the following relations between observables 
0 ^( 02 ) and the helicity coefficients: 


0^^\a2) = 

0 f)(02)= ( 


\da2' 




Vdo 2 / 


d0 


de^ 



hie, ] 

- -\a„ 

Re(7ljp7ll) 

da2d9i 

J da2d9i) 
1 

Epl^),P 

, d^T 

fee- 1 

3 . 

1 A n 

Re(7ljp7l'j) 

^ da2d92 

J ^^^da2d92) 

1 2 



(26) 


under the condition G ^ 0 ,mx — yjnri^z + ^2 • 

One can write these two formulas in the following way: 




J da2d9i^2 J da2d9i^2' 2 


Re(Glj|Ml) 


(27) 


Then we deduce that 




r+/3 


d 6»2 


da 2 d 62 y ^ ^ ^ \da2^ 

^+Q: 


J d92...) = (^) '( / d9. 


d^T 

da2d9i 


\+P 


d9^ 


f+a 


= - ^(sin,5 — sin a) (3 + sin^ a + sin^ (3 + sin a sin/3) 

lyl'P \ 

+ 3 ° (sin,8cos^/3 —sinacos^a)). 


'Epl^p 
0<a</3< |, 


(28) 


2 2 




dcL2d0\d02 


da2d6 




0 0 

- - A J. An II _ 

16 Epl^),P ’ 


(29) 


dr\-i 




dip 


d^T 

da2dp 


— / dp... + / dp... — / dp... + 


dp..)j 


1 

2^ Epi^;,P 


(30) 
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2 


d^T 

da2d^p 


— / dif...+ / dip... — / dip...] = — 


1 Im(AjpAl) 


(31) 


2 


da2dip 


2'jr 


— / dip... + / 


_ 9 ^ ^ Re(A'Mj|) 

32 /i /. ^ 


(32) 


TT 2 71' ■^TT 


= qV^^/ 2,1 


0 0 

lm(A'*Aj|) 


(33) 


O 




0 

3 z:^, Re{A'Q^A'_^) 

Epl^(,P ’ 


d^r 


(^J d,ip(^J d 1.2 J ^ da2d9i^2dip 

c 

= 


2 '• ZTT 

) - / + J dip..^ 


27r 




Im(AyAl) 
Epl^(,P ■ 


(34) 


(35) 


From the measured observables 0 ^( 02 ) one can get constraints on the dependences of the couplings 0 ^( 02 ), h'^{a 2 ) 
and € 2 ( 02 ). As for 02 ( 02)5 it can be measured at a fixed value of 1^/02 and at various values of the parameters 
/3 and a. Then, after obtaining central values and uncertainties of a quantity |Agp/^p|App from Eq. (1^!^ at 
several sets of values of /3, a, one can combine these central values and uncertainties and thereby get a value of 
l^oP/Sp l^pP with greater precision than in case of any particular values of /3, a. 

As an illustration of the behavior of these observables, in Fig. 0] we show their ^^-dependence with the constant 
a'hz, b'f^z’ ^'hz from the sets (IMll . The observable 02 ( 02 ) is presented for (3 = 90° and a = 70°. 

As one can see, for each of the observables 02 ( 02 ), 03 ( 02 ), 04 ( 02 ), 05 ( 02 ) their dependences on y/a 2 for all the 
four sets (EH) are very close. The observables 02 ( 02 ) (at /3 = 90°, a = 70°) and 04 ( 02 ) are relatively large with 
the maximum values greater than 0.05, and thus to measure these observables a relatively small amount of data is 
needed, while 03 ( 02 ) and 05 ( 02 ) are smaller, which complicates their experimental observation. 

Further, Oj^’^^(o 2 ), 05 ( 02 ), Og^’^^(o 2 ), 09 ( 02 ) vanish for €'^ 2 ( 02 ) = 0, according to Eqs. (ITfl) . (I^Tl) . (1541) . (1551) 
and Fig.m Therefore, these observables can give significant constraints on the CP-odd coupling €(^ 2 ( 02 ), although 
their moduli are relatively small. 

The functions 07 ^’^^(o 2 ) are proportional to lm.{a'^zb'hz) (®®® (l33l) . ([ 6 ])), and, consequently, they are equal to 
zero for any set from (1241) . Among all the observables under consideration, 07 ^’^^(a 2 ) are the only ones vanishing 
in case b’^zi^^) = 0 for any 0 ) 12 ( 02 ) and c'f^zi'^'^)- Therefore, knowing the dependences ’ { 0 , 2 ) allows one to get 
notable constraints on the function 65 ^ 2 ( 02 ). Although in case (15511 these observables turn out to be relatively small 
in absolute value (see Eig. [5D. 
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VoT, GeV 
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Fig. 4: The observables 0^'^\ O 2 (at /3 = 90°, a = 70°), O3, O4, O5, Oe, O9 for the decay h —>■ —>■ 

^ 2^2 ih = e)A‘) 2 ")^i / ^ 2 ) as functions of in case |a)iz| = = 0 , c'/,^ = 0 (solid lines); a'f^z = 1 , 6 (,z = 0 , = 

0.5 (dashed lines); a'f^z = l,5!iz = 0, = 0.5i (dash-dotted lines); a(,z = l,ti)iz = ~0.5, = 0 (dotted lines). 



V^, GeV 


Fig. 5: The observables for the decay h —>■ > /j lfl 2 it ih = 65 /r, r, / I 2 ) versus y/^ in case a'f^z = 1 ) ^Iz — 

—0.5i, c'f^z = 0. 
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Note that from (1771) - (155)) . regardless of the values of the couplings 02 ( 02 ), ^ 2 ( 02 ) and € 2 ( 02 ), it follows that for 
any 02 


e 


r 3 , 

[ 4"^/i,2 


r^/1,2], O2 G [ 0 , 1 ], O3 G [ 0 , —Af^Af^], 04,05 G[- 


— —] 
27r’ 27r^’ 


06,09 G Av/27r%%], 0('’^\0^^) G (36) 

Since A^.- = A^~ = A^- « 0.214, = A^^ = = 1, = Ac = At « 0.697, Ad = Ag = At, k. 0.941, the 

moduli of 02 ^’^^(a 2 ), 07 ^’^^(a 2 ), Og^’^^(a 2 ), 03 ( 02 ), 05 ( 02 ), 09 ( 02 ) for the decays (P) with quarks and/or neutrinos 
in the final states are greater than those for the decays o to leptons, therefore, the former processes seem more 
feasible for experimental study. On the other hand, detection of leptons is much simpler. That is why the study 
of each decay channel of the type o has advantages and disadvantages which strongly depend on experimental 
methods and parameters of detectors. Consequently, measurement of the observables 0\ ’ ( 02 ), ..., 09 ( 02 ) for 
various decay channels and for various invariant masses of the fermion pair {^/a^) may help to put constraints on 
the XZZ couplings 02 ( 02 ), ^ 2 (^ 2 ) and 02 ( 02 ). 


III. CONCLUSIONS 


In the present paper the decay of a neutral particle X with zero spin and arbitrary CP parity into two off-mass- 
shell Z bosons {Zi and Z|) each of which decays to a fermion-antifermion pair, i.e. the decay X Z\Z 2 —>■ 
/ 1 / 1 / 2 / 2 , has been considered. The given decay has been examined at tree level for the non-identical fermions, 
fi ^ f 2 - In the approximation of the massless fermions a formula for the fully differential width has been obtained. 
It has been established that the narrow-Z-width approximation is applicable for finding differential decay widths 
oi X ^ Z*Z 2 ^ / 1 / 1 / 2/2 only if the invariant mass of the pair f 2 f 2 lies in an interval ^ 0 , nix — \/ + ^2 ■ 
If the parameter £2 gets larger, the accuracy of the used approximation increases, but the interval in which the 
approximation is valid reduces. As an optimal value of £2 we have chosen £2 = Sm^Tz- 

In the narrow-Z-width approximation, but without the neglect of T^ in the propagator of Z|, a formula for 
the total width of the decay o and the total width of /i —>■ Z/Zf have been derived. The former formula is 
valid in case /i = /2 as well. Note that in Ref. 23| within the framework of the SM the total width of the decay 


X —>■ ZZ* —>■ Zjj has been found in the approximation Vz ~ 0 in the propagator of Z*. In an analogous way one 
can obtain the total width of the decay (HD in the SM after the neglect of T^ in the propagator of Z 2 , however the 
formula ()A2|) . derived in the present paper, is more general and more precise. 

Using the CMS data js), we have found constraints on the couplings a/ 2 > which determine the hZZ 

interaction and the CP properties of the boson h detected in the experiments l|. Comparing our constraints with 


those reported in Refs. 


36 


37|, one can notice appreciable overlaps between the three results. Besides, we have 


derived an allowed interval for a ratio not studied in 361137|. Taking our allowed regions into account, we have 
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selected several sets of values of the couplings (<Z = c) and analyzed results for these sets. 

The observables 0\ ’ ' ( 02 ), 09 (a 2 ), measurement of which will allow one to get constraints on the dependences 

of q'^ on y/a 2 , are defined. It is shown that the observables 05 ( 02 ), Og^’^^(a 2 ), 09 ( 02 ) become zero in 

case 0 ^( 02 ) = 0, and therefore their experimental dependences on can put significant constraints on the 
CP-odd coupling 0 ^( 02 ). The observables 07 ^’^^(o 2 ) vanish if b'^{a 2 ) = 0, and, therefore, their measurement is 
important for finding the OP-even coupling b'^(a 2 ). 

Note that the absolute values of 0^i'^\a2), 0^'^\a2), 0^'^\a2), 03 ( 02 ), 05 ( 02 ) and 09 ( 02 ) for the decays ([T]) 
where fi and/or /2 is a quark or a neutrino are greater than those for the processes in which the fermions are 
leptons. At the same time, the processes with the leptons are much more convenient from the experimental point 
of view. 

Thus, measurement of the observables o[^’^\a 2 ), ..., 09 ( 02 ) for the decays ([1]) can help to clarify the CP 
properties of the particle X and the structure of the amplitude of the decay X —>■ Z/ . 

The authors thank Sergiy Ivashyn for useful discussions. The work is partially supported by the National 
Academy of Sciences of Ukraine (project 11,0-15-1/2015) and the Ministry of Education and Science of Ukraine 
(project 0115U00473). 

Appendix A: Calculation of the total widths of the decays X ^ Zi Z2 ^ / 1 / 1 / 2/2 and h ^ ZIZ2 

In this Appendix we calculate the total width of the decay X ^ ^ / 1 / 1 / 2/2 for the my-values such that 

mx > rnz and for the dependences of 0 ^( 01 , 02 ), 5z(ai,a2), 0 ^( 01 , 02 ) such that the differential width has 

a sharp maximum when = mz or ,^02 = mz and the functions jczl, cos(arg5'2 — arga'^) (q'z = 

02 ); q = a,b,c) are independent of 02 . For example, \az\ ^ 1, bz ^ 0, Cz ^ 0 are such dependences (see 
Sec. EH). Then we calculate the total decay width oi h ^ Z^ Z 2 and examine the applicability of an approximation 
ft! 0 for derivation of the total widths. 

1 . The total width of the decay X ^ Z*Z2 / 1 / 1 / 2/2 

Analogously to the derivation of Eq. (11111 . we find that 

275 f 

Tft-— / da2/(m|,a2). (Al) 

mzFz J 
0 

Let us consider the case wherein \a'^\, \b'^\, cos(arg6'^ — arga^) are independent of 02 . Having exactly 

calculated the integral in Eq. (lAll) with allowance for Eqs. ciiid (|6j) . wG obtciiii! 


r ft!(aj^ -I- fo '^/ 2 )/o(®z, Cz, ^z,^z, s), 


(A2) 
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where 


/o(azi s) = 


y/2G%.m7zmx 


2i0337r4r2 

(1 - a) (-24(23a - b)\a'z\^ + (Za^ - 37a^ - a(235 + 6/3^) + 77 - 54:P^)\b'z\‘^ 

- 16(2a2 + 26a - 13 + 3^^)Re{azb'z) + 64a{a^ + 40a - 11 + QP^)\c'z\'^^ + 61n 

X (4(12a^ - 18a + 3 - I3^)\az\'^ + (30a2 - 10a(3 - /3^) + 5 - 10^^ + I3‘^)\b'z\'^ 

+ 8(6a2 - a(9 - /3^) + 2 - 2l3^)Re{azb'z) - 32a(6a2 - a(9 - /S^) + 1 - 3l3‘^)\c'z\^^ 


3^2 , 


P(a, /3, a'^, 6^, c^, r+, -4/3r_) 


X In 


(1 — a)^y/(4a — 1 + /3^)^ + 4/3^ + (3a — 1)^ + /3^(a + 1)^ + S'/2{1 — a)((3a — l)r_ — /3(a + l)r+) 


4a (a^ + /3^) 


+ 2P(a, /3, a'z, b'z, c'z, r_, 4/3r+) - arg(-a(3a - 1 + /3^) - /3^ + s-^(/3r+ - ar_) 


+ i(l - - P{1 - a) 


V2 




— mzT z 

- -77^’ 


P{a, /?, a'z, b'z, c'z, x, y) = 2(2a;(12a^ — 4a + 1 — 0^) + y(6a — l))|a 2 p + (a;(16a^ — 8a(l — 0^) + 1 

- 6/3^ + P0 + j/(4a - 1 + P0)\b'z\'^ + (4a;(8a^ - 2a(3 - /3^) + 1 - 3/3^) 
+ y(8a — 3 + /3^))Re(a262) ~ 8a(4x(4a^ — a(l — 0^) — 2fi0 

+ y(6a-l + /3^))l4P, 


(A3) 


(A4) 


r± = 


: - 1 + /32)2 + 4/32 ± (4a - 1 + /32). 


(A5) 


In place of s one may take 1 or -1 {fQ(a'z,b'z,c'z,Tnz,7^z,—^) = fo(p-'z,b'z,c!z,rnz,7^z,7)). In this article the 
argument arg 2 of a complex number 2 is defined as follows: 


Im z 

arg 2 = arctan-h 7rn(Re 2 , Im 2 ) V 2 £ ClRe 2 ^ 0, 

Re 2 

arg 2 = TT + 0(—Im 2 )^ V 2 £C|(Re 2 = 0 and Im2^0), 


(A6) 


where n(x, y) = Q{—x) + 2 0(x)0(—y) Vx 0 0, 


0(x) = 0 VxG(— oo,0], 0(x) = 1 VxG(0,+oo). (A7) 

From the definition (IA6I) it follows that arg 2 is the angle counted clockwise on the complex plane from the vector 
(Re 2 , Im 2 ) towards the vector (1, 0) and arg 2 £ [0, 27r). Sometimes in literature a different function 


arg' 2 = arg 2 — 27r0(—Im 2 ) 


(A8) 
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is used as the argument of z. From (IA 8 I) it follows that arg'z € (— 7 r, 7 r]. Note that we have already used argz 
above in the expression cos(arg 62 — argaz), but since cos(arg 6 z — aigaz) = cos(arg' 6 z — arg'a^), at that point 
the distinction between argz and arg' 2 was irrelevant. 

Calcnlating the integral over 02 in Eq. jsa, one hnds an antiderivative of f(rri^,a 2 ) on the interval 
[0, {mx — In this antiderivative the fnnction argui(a 2 ) natnrally appears, where Mi(a 2 ) is a complex-valued 

dimensionless function snch that 


Va2 G [ 0 , (mjf - mz)^) Im 141(02) 7^ 0 , 

iTnui^irax — mz)^) = 0, Reui((mx — mz)^) < 0. (A9) 


arg' 01 ( 02 ) does not emerge in place of argiti(a 2 ) since, according to (IA8I) . the function arg'z has a discontinnity 
on the half-line Imz = 0,Rez < 0 and thus arg' 01 ( 02 ) has a discontinuity at the point 02 = {mx — mz)^■ To 
avoid this drawback it is convenient to use arg z in Eq. (IA3I) . 

Note that in case of the Higgs boson, i.e. X = h, Eq. (EH) can also be written in terms of the fnnction arg': for 
this one has to substitute in Eq. (IA3I) tt — arg... by stt — arg'..., since according to (IA8I) and to data of Table U 
TT — arg ... = STT — arg'.... 

In case of the identical fermions, /i = / 2 , one may neglect the interference term and then in order to obtain a 
formnla for E one has to mnltiply the right-hand side of the relation (IA2I1 by (hi view of the identity of the 
final fermions) and by 2 (since the contribution of the diagram with the permutation of the particles to T is equal 
to that of the diagram without the permutation), i.e. to multiply the right-hand side by Consequently, for any 
fi and /2 

T Ri(l — -< 5 /i/ 2 )(a/i + i'/i)(®/2 + ^/2)/o(oz: b'z, c'2, mz, Tz, s) = Tr^, (AlO) 


where = 0(1) at /i 7 ^ /2 (/i = / 2 ). The neglected interference term seems small based on qualitative 


arguments of Ref. 


38| . For a quantitative estimate we can nse Ref. 


30| (see Table 1 there), according to which the 


interference contribution to T{h Z^Z 2 —>■ 4e) in the SM at tree level is 5.80% for mt = 140 GeV. 


In Ref. 


23| the width of the decay h —>■ ZZ* —>■ Zff has been derived at tree level in the SM after the neglect 


of in the propagator of Z*. Following 


23l |. when calcnlating the integral in Eq. (lAll) . in the expression for 


f{m^,a 2 ) we may also neglect F^, and then we obtain the following approximate formnla for F in the SM: 
r|sM ~(1- 2^A/2)——(a/i +^/i)(“/2 +^/2) 


, „ 1 — 8a -I- 20 a^ 

X 6- , — arccos 

V4a- 1 


3a- 1 

2ai 


-——{2 — 13a -I- 47a^) -I- 3(1 — 6 a -I- 4a^) In — 
a a 


= r, 


OlSM- 


(All) 


From (lAllI) and (IA10|) we obtain that at mx = mu 


TqIsm ~ 1.001 X Ty-z\sm- 


(A 12 ) 
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Besides, Fq > Fr^ (for any /i,/2) since when deriving the formula for Fq one neglects the width Fz 

in /(m|,a2) and the value of the integral increases. Still according to (IA12|1 . the difference between FqIsm and 
BrzIsM is about one per mille. 

Finally, note that at mx = rrih we can represent the dependence of the function /o on the XZZ couplings 
a'z,b'z,c'z in the convenient form: 

foia'z,b'z,c'z,mz,Fz,s) Ki (3.359|a2p + 0.05216^^ + 0.594Re(a'^6'2) + 0.125|c2p) keV. (A13) 


2 . The total width of the decay h ^ ZIZ 2 


The total decay width F(/i —>■ Z^Z^) is 

F(h ^ ziz;) = EE ^Imx^rrih ; (A14) 

h h>fi 

where the sums run over the fermions e“, r“, r'e, Ui, Ci, di, Si, bi (since mh S (4m{,,2mt]), i = r,g,b is 

an index of quark color. It follows from Eqs. IIA14|1 . (IA12p that in the SM 


Fo(h ^ Z*Z*) ft! 1.001 X Fr^(h ^ ^ 1 *^ 2 )- 


(A15) 


Further we use Eq. (lAlOl) since it is more precise than Eq. (lAllll and consider the case wherein \a'f^z\j Whzl’ Whz\^ 
cos(arg5'^2 ~ argaj^^) do not depend on 02. From Eqs. (IA14|) . (lAlOl) . (|A3p we derive that 


T{h —>• ^*^2 ) ~ /o(«/iZj mz, Fz, s)|mx=m^ ( 2 E^“/i 

V fi 


1 


fi 
/o(---) 


9 E E E(“/+^/) “ 


h / 2//1 

/o(-) (103 
18 \ 2 


- 100 


/ mw 
V mz 


80 


mw V 

mz ) 


Carrying out calculations, we find the total decay width for the sets (IMl) . (P5|) : 


(A16) 


Whz\ = 1, Kz = 0, c'^z =0^F{h^ ZIZ*) ft 91.16^11;|6 keV, 

Whz\ = 1, Kz = 0, \c'hz\ = 0-5 ^ F(h ^ ZlZl) ft 92.0i;il«? keV, 
a'hz = 1, Kz = -0.5, c'^z = 0 ^ F(h ^ ZIZ*) ft 83.45^15;?® keV, 

a'hz = 1, Kz = ±0.5z, c'^z = Q^F{h^ Z^Z*) ft keV. (A17) 


The uncertainties shown in Eqs. (|A17I) are calculated by finding the maximum and minimum values of the 
function F{h —>■ Z^Z^) in the region v € [uq — 3(T„,no + 3cr„] (v = GF,mh,mz,mw,Fz)- Here vq is the central 
value of a quantity v, tT„ is the 1-standard-deviation uncertainty of v; according to the data of Table III Gfo = 
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1.1663787 X 10-5 GeV’^, crcp = 6 x lO’^^ GeV^, ruho = 125.7 GeV, = 0.4 GeV etc. 
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